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Abstract 

Riemann and sectional curvatures of magnetic twisted flux tubes in Riemannian manifold are 
computed to investigate the stability of the plasma astrophysical tubes. The geodesic equations 
are used to show that in the case of thick magnetic tubes, the curvature of planar (Frenet 
torsion-free) tubes have the effect ct of damping the flow speed along the tube. Stability of 
geodesic flows in the Riemannian twisted thin tubes (almost filaments), against constant radial 
perturbations is investigated by using the method of negative sectional curvature for unstable 
flows. No special form of the flow like Beltrami flows is admitted, and the proof is general 
for the case of thin magnetic flux tubes. In the magnetic equilibrium state, the twist of the 
tube is shown to display also a damping effect on the toroidal velocity of the plasma flow. It 
is found that for positive perturbations and angular speed of the flow, instability is achieved , 
since the sectional Ricci curvature of the magnetic twisted tube metric is negative. Solar flare 
production may appear from these geometrical instabilities of the twisted solar loops. PACS 

numbers: 02.40.Hw:Riemannian geometries 
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I Introduction 



The stability of geodesic flows have been recently investigated by Kambe [1] by making use of 
the technique of Ricci sectional curvature [2], where the negative sectional curvature indicates 
instability of the flow, while positivity or nul indicates stability. In the case of instability the 
geodesies deviate from the perturbation of the fluid. Following the work of D. Anosov [3] on 
the perturbation in geodesic flows in three-dimensional Riemannian geometry, in this paper 
the sectional Riemann curvature of the geodesic flow for a Riemannian flux tube [4, 5], where 
the axis of the tube flow possesses Frenet curvature and torsion. In the approximation of a 
thin tube where the radius of the tube is almost null, we show that the flows are unstable, 
against orthogonal perturbations, which is equivalently due to the negativity of the sectional 
Ricci sectional curvature. Throughout the paper the ellegant coordinate-free language of 
differential geometry [6] is used. The paper is organized as follows: Section II presents a 
brief review of Riemannian geometry in the coordinate free language. Section III presents the 
geodesic flow computation of the Christoffel symbols for the thick flux tube, where we show 
that the curvature of the tube axihe speed of the flow. Section IV presents the computation 
of the instability of Riemannian tube flow. Section V presents the conclusions. 

II Ricci and sectional Riemann curvatures 

In this section we make a brief review of the differential geometry of surfaces in coordinate- free 
language. The Riemann curvature is defined by 



where XeTAi is the vector representation which is defined on the tangent space TAi to the 
manifold M.. Here Vx^ represents the covariant derivative given by 



which for the physicists is intuitive, since we are saying that we are performing derivative along 
the X direction. The expression [-^,1^] represents the commutator, which on a vector basis 



R{X, Y)Z := VxVyZ - VyVxZ - V^x,y]Z 



{111) 



VxY = {X.V)Y 



(11.2) 
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frame e; in this tangent sub-manifold defined by 

X = Xuek (II.3) 

or in the dual basis 

X = X^dk (II.4) 

can be expressed as 

[X,Y] = {XXfdk (11.5) 
In this same coordinate basis now we are able to write the curvature expression (II. 1) as 

R{X, Y)Z [R'^upZ'X''YP\di (II.6) 

where the Einstein summation convention of tensor calculus is used. The expression Y)Y 
which we shall compute bellow is called Ricci curvature. The sectional curvature which is very 
useful in future computations is defined by 

where S{X, Y) is defined by 

s{x,Y) := < x,r >2 (ii.s) 

where the symbol <, > implies internal product. 
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Ill Geodesic equations in Riemannian tube metric 

In this section we shall consider the twisted flux tube Riemann metric. The metric g{X,Y) 
line element can be defined as [4, 5] 

ds"^ = dr^ + r^deR^ + K\s)ds'^ (III.9) 

This line element was used previously by Ricca [4] and the author [5] as a magnetic flux tubes 
with applications in solar and plasma astrophysics. This is a Riemannian line element 

ds^ = Qijdx'dx^ (III. 10) 

if the tube coordinates are (r, Or, s) [4] where 9{s) — Or — f rds where r is the Frenet torsion 
of the tube axis and K(s) is given by 

K^{s) = [1 - rK(s)cos^(s)]2 (III.ll) 

Let us now compute the geodesic equations 

— + rjkv'v' = (III.12) 

where — ^ and — ^ the Riemann-Christoffel symbols are given by 

rV = + 9ik,j - 9jk,i] (111.13) 

The only nonvanishing components of the Christoffel symbols or Levi-Civita [6] connections 
of the flux tube are 

= - (111.14) 

r 

r^33 = z^^i^l^ (III.15) 

r 

r^i = ^ [Tsm9 + -cos9] (III. 16) 

2K K 

r\s^K-'K' ^-r\, (III. 17) 

A simple example can be given by writing the geodesic equation for the untwisted tube, where 
vo — 0, as 

s + r^ip-s^] = (III.18) 



where substitution of the component F^n above yields 

dlnvs = -d{lnK) (111.19) 
where we use the chain rule of differential calculus, ~ — . Solution of equation (111.19) 

etc Vs 

is 

Vs = vqk-^ (III.20) 

where vq is an integration constant. This solution tells us that when the curvature tends to 
oo the velocity along the tube axis vanishes. Physically this means that the curvature acts as 
a damping to the flow, along the tube. In the next section we investigate in some detail the 
stability of an incompressible or volume preserving flow, using the method of the sign of the 
Ricci sectional curvature. 



IV Sectional curvature and plasma flow stability 

One of the most important features of the investigation of the stability of flows in the Euclidean 
manifold E^, is the comprehension of the fact that the covariant derivative in the flow curved 
manifold is given by the gradient operator in curvilinear coordinates. Thus to compute the 
Riemann sectional curvature above, we need to make use of the grad operator in the twisted 
flux tube Riemannian metric given by 

V =[dr,r-^dg^,K-^d,] (IV.21) 

Since the axis of the tube undergoes torsion and curvature, we need some dynamical relations 
from vector analysis and differential geometry of curves [7] such as the Frenet frame {t,n,b) 
equations 

F = (IV.22) 

n' = -Kt + Tb (IV.23) 

b' = -rn (IV.24) 

and the other frame vectors are 

= ncos9 + bsinO (IV.25) 
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eg = —fisinO + bcosO (IV. 26) 

deee = -n[{l + T'^K)sine + cosO] - b[cose + sinO] (IV.27) 
Let the constant perturbation be given by 

X = Ur^er (IV.28) 

The upper index one in this expression refers to the fact that the background original value of 
Ur was considered as = to form the tube. The other variable Y is given by 

Y = ueee + u^t (IV.29) 

Therefore to compute the Ricci tensor step by step we start by the term 

VxY ^Ur^^'^dr[ueee + Ust\ (IV.30) 

which vanishes since we addopt here the approximation Ur^^^drlue] ~ along with the same 
relation to the radial partial derivative of Ur. So 

VyVxY ^ (IV.31) 

Now the second term in the Ricci tensor is 

Vx^yY = -Ur^^\-\g[^^-^{ncose + bsinO)] (IV.32) 

T 

where we have used the approximation of the thin tube where K{s) ~ 1 and r ~ 

[X.Y] = Ur^^^ [r-^ug - TUs] [e^ - T-H"! (IV.33) 

which implies that 

V[x,y]l^ = Ur'^^\r-\g - TUs][-T-\gdsee + nn] (IV.34) 

The Ricci tensor is 

R{X, Y)Y = -Ur^^^-^ugT-^dsee (IV.35) 
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In the previous computations we have made use of the imcompressibility of the flow 

V.u = (IV.36) 

which is 

dgUe = ttkuq fti (IV. 37) 

since r on the RHS of equation (IV.37). The sectional curvature is thus 

when the tube is strongly twisted, uq^ » thus the sectional Ricci curvature is 

<R{X,Y)Y,X> \\ + t{s)kcosB\ 
^ S{XX) ^ ru^'\ue ^^^"^^^ 

when the tube, besides is planar or torsion vanishes the last expression reduces to 

Note that when both angular velocity and perturbation both keep the same sign, the sectional 
curvature K{X^ Y) is negative and the flow along the Riemannian flux tube is unstable. There 
is singularity in this sectional Riemannian curvature in r ~ 0. Note that if we consider the 
imcompressibility equation (IV.36) as 

—ue = ue—sin9 (IV.41) 

Since the magnetic field B is also divergence- free, where B,. vanishes in the equilibrium state, 
the same equation for the poloidal magnetic field component Bq is obtained 

—Be = Be'^sinO (IV.42) 
OS K 

Putting K ^ for the thin tube approximation, and the equilibrium state magnetohydrody- 
namics (MHD) equation 

Vx[mxS]=0 (IV.43) 

and taken the most simple solution 

uxB^O (IV.44) 
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one obtains 

ugB, = UsBe (IV.45) 

since Ur and Br both vanishes before perturbation in the radial direction of the magnetic flux 
tube. Taken into account the relation derived by Ricca [4] on the ratio between the poloidal 
and toridal components of the magnetic field 

where we have taken K = 1. Substitution of this expression into (IV.45) and using the 
equation for the divergence-free of the velocity plasma flow field yields 

^ m -cose] 

27rrTw ^ ' 

where 

2'KrTwUs /„ . , X 

ue = (IV.48) 

or ue — 1 — cos9\ which upon substitution in the sectional curvature (IV.40) yields 

which shows that since cosOpel, the sectional curvature is singular or negative which implies 
instability of the solar loops, as long as the radial perturbation on the twisted solar loop is 
positive. In other words, if the tube is radially expanding as happens on the suface of the sun, 
the solar loop is highly unstable, which is reasonable for the production of solar flares. 



V Conclusions 

An important issue in plasma astrophysics as well as in fluid mechanics is to know when a fluid, 
charged or not, is unstable or not. In this paper we discuss and present incompressible flows 
and investigate their stability. Instability is obtained even before the singularity is achieved. 
Inflexional desiquilibrium of solar loops has been also recently investigated by Ricca [4], based 
on the above Riemann twisted magnetic flux tube. He proved a theorem where inflexion 
desiquilibrium is proved from a state of MHD equilibrium during passage throghout inflexional 
points in the solar loop. In this sense, though obtained from completely different methods it 
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seems that coincide with those of Ricca's [4]. Solar flare production may appear from these 
geometrical instabilities of the twisted solar loops. 
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